Towards unified understanding of conductance of stretched monatomic contacts 
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When monatomic contacts are stretched, their conductance behaves in qualitatively different 
ways depending on their constituent atomic elements. Under a single assumption of resonance 
formation, we show that various conductance behavior can be understood in a unified way in terms 
of the response of the resonance to stretching. This analysis clarifies the crucial roles played by the 
number of valence electrons, charge neutrality, and orbital shapes. 

PACS numbers: 73.40.Cg, 73.40.Jn, 73.63. Rt 
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Monatomic contacts with one-atom-wide cross-section 
(referred to as "contacts" hereafter) have been realized by 
scanning tunneling microscopes or mechanically control- 
lable break junctions 0. Upon stretching, the conduc- 
tance G of various contacts shows different behavior. In 
contacts made of monovalent atoms such as Na and Au, 
only one channel contributes to transport and G stays 
at the quantized value G = 2e 2 /h during the stretch- 
ing. 0, 03, [J IE El In contacts made of polyvalent atoms 
such as Al and Pb, on the other hand, current is car- 
ried by three channels and G does not stay at quantized 
values during the stretching 0, El ED ■ Thus monova- 
lent and polyvalent contacts are qualitatively different. 
Moreover there are differences even between polyvalent 
contacts; during the stretching, G increases (decreases) 
for Al (Pb) contacts. 

The different conductance features were addressed by 
numerical calculations 0, 0, 0, 0, Q, 0] and some 
calculations reproduced experimental results for certain 
contacts. However the calculational results are often too 
specific to the particular contacts under calculation and 
the understanding on the origin of the aforementioned 
difference is yet far from satisfactory. For example, the 
connection between the conductance features and the 
number of valence electrons is poorly understood. In 
this respect, it is desirable to have an alternative ap- 
proach that allows easier comparison of various contacts 
and thus provides an insightful explanation of the dif- 
ference. We also remark that numerical calculations of 
conductance are sometimes subject to strong finite size 
effects [l^ . Thus a naive comparison may be dangerous 
between experimental findings and numerical calculation 
results based on small cluster modelings of atomic con- 
tacts. For the behavior of Al contacts, existing calcu- 
latio ns IToL Il2j have reported qualitatively different re- 
sults [lTj and we suspect that the difference may be re- 
lated to the finite size effect. 



In this paper, we present an analysis that allows one to 
capture the connection between the conductance prop- 
erties and microscopic features such as the number of 
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FIG. 1: A model of a monatomic Al atom contact |T3 |. 



valence electrons, charge neutrality, and orbital shapes. 
This connection provides a unified understanding of var- 
ious conductance behavior. This analysis is applicable, 
provided that electron transport has resonant character. 

The assumption of the resonance is very plausible in a 
stretched contact. First of all, the central atom of a con- 
tact (Fig.ni El 

is relatively weakly coupled to its envi- 
ronment even without stretching since it has a fewer num- 
ber of nearest neighbors. When stretched, the coupling 
becomes weaker and it is natural to expect that atomic 
orbitals of the central atom give rise to the resonance EH • 
Among the three contacts (Na, Pb, Al) on which we are 
going to concentrate, the resonance formation is alr eady 
demonstrated for Na and Pb contacts [ToL Iu3 . Hil Il5| . 
For Al contacts, we present in this paper ah initio calcu- 
lations that confirm the resonance formation. 

Transport channels.- The number of transport chan- 
nels depends on the number of valence electrons 0, El ; 
while monovalent contacts have only one transport chan- 
nel, polyvalent contacts have multiple channels. We be- 
gin our analysis by identifying "eigen" -transport chan- 
nels, which are "eigenvectors" of a proper transmission 
matrix and decoupled from each other in this sense. The 
identification is straightforward for a contact with a sym- 
metric shape; recent experiments [2(| showed that con- 
tacts tend to align their axes with high symmetry lattice 
directions so that contacts have a symmetric shape. As 
an illustration, we describe the identification of 3 chan- 
nels E3 f° r a (HI) Al contact with the 2n/3 rotation 
symmetry and the inversion symmetry (Fig. QJ. The 
rotation symmetry with respect to the z-axis (the hor- 
izontal axis in Fig. ^) makes the angular momentum m 
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(= 0, ±1) under the 2ir/3 rotation a good quantum num- 
ber. Then a traveling wave packet with a given m can 
be scattered into states with the same m only and thus 
the three eigen-transport channels of Al contacts can be 
characterized by m. There is a natural link between the 
channels and the orbitals of the central Al atom; m = 
channel couples with 3s and 3p z orbitals, while m = ±1 
channels with 3p x ±i3p y . Here the factor ±z implies that 
the two channels m = ±1 are mutually related by the 
time reversal and thus degenerate. 

Friedel sum rule.- The conductance of a contact is 
given by G = G J2 m T m(E F ), where T m {E F ) is the 
transmission probability of the channel m at the Fermi 
energy Ep. In order to examine the connection between 
T m and microscopic features, it is useful to use the Friedel 
sum rule [2ll 12*2^ . For a contact with the inversion sym- 
metry, which makes the parity under the inversion a good 
quantum number, the Friedel sum rule [2l| results in 



Tm{Ep) = sin 2 ( — AN, 



(1) 
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where AN m = 

2j Ep dEp m ^ )(E), and p m . e / is the density of 
states (DOS) with even/odd parity in the channel m 
(spin degeneracy assumed). 

Resonant transport and charge neutrality. - We then ex- 
amine how AN m 's are affected by microscopic features. 
For the resonant transport, AN m is governed by reso- 
nances, AN m = 2(R m , e — R m ,o), where R m ,e(o) is the 
number of even (odd) parity resonances in the channel m 
below Ep (2 from spin degeneracy). Partially filled reso- 
nances (near Ep) give fractional contribution to R m .e/o- 

Another important relation comes from the charge 
neutrality, which is an excellent approximation in metal- 
lic systems 0, ^| . The charge neutrality near the cen- 
tral atom results in 2 ^ m (i? m , c +fl m ,„) = N', where N' is 
the number of valence electrons in the central atom [2^ . 

To make a further progress, we first note that T m is 
periodic in R m ,ejo with the period 1 [Eq. JJJ]. We thus 
ignore the integer part and focus on the fractional part 
of R m ,e/o arising from the partially filled resonance. For 
contacts under consideration, moreover, central atoms 
possess only a few atomic orbitals and as a result, there 
is at best one partially filled resonance in each chan- 
nel. This is verified by our ab initio calculation given 
below (Al contact) and by Refs. [13,113 (Pb and Na con- 
tacts, respectively). Then for each m, either R m , e — or 
Rm,o = 0, and the new equality \AN m \ = 2(R m , e + R m ,o) 
holds. The charge neutrality condition then becomes 



|AiV m |=^ ff , 



(2) 



where N' cS can differ from N' by an integer multiple of 2 
due to the neglect of the integer part of R m . e jo- 
Mono, vs. polyvalent contacts.- For monovalent con- 
tacts such as Na [l4|, where there is only one transport 



channel and N' = N^ B = 1, Eq. © reduces to |AiV| = 1. 
According to Eq. (Q), Eq. (J2J then implies G — G re- 
gardless of other details such as the degree of stretching, 
explaining the quantization of G in experiments. 

In polyvalent contacts, on the other hand, there are 
multiple channels and multiple AN m 's. Then the single 
constraint alone cannot completely determine AN m 's 
and G needs not be quantized, explaining the difference 
between monovalent and polyvalent contacts. This is one 
of the main results of this paper. 

As a concrete example of polyvalent contacts, we ex- 
amine the symmetric Al contact. When its central atom 
is weakly coupled to the electrodes, its atomic orbitals 
3s, 3p x , 3p y , 3>p z give rise to one (almost) fully occupied 
resonance in the m = channel (mostly s character) 
and three partially filled resonances, one in each chan- 
nel. This occupation of resonances can be easily inferred 
from the knowledge of an isolated Al atom. When the 
two electrons in the fully occupied resonance are ignored, 
N' eS is N' — 2 = 1 and from Eq. @, we have 
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■2|AJV±i| = 1, 



(3) 



where the degeneracy of the m = ±1 channels is used. 
From Eq. JJ), the conductance then reads 



cos (irAN ±1 ) 



(4) 



Thus G varies as AN±± changes during the stretching. 

Orbital shapes.- The value of AiV±i is influenced by 
the shapes of s and p orbitals, which make the ppix cou- 
pling (responsible for the partially filled m = ±1 reso- 
nances) weaker than the ppo~ or spa coupling (respon- 
sible for the partially filled m — resonance). Con- 
sequently |AJV±i| < |A7V |, and from Eq. 10), we find 
|AiV±i| < 1/3. The stretc hing magnifies the effect of 
the orbital shape difference [lfj, |24j and thus decreases 
|AAT±i| further below 1/3. Then according to Eq. (J2J, G 
increases monotonically durin g th e stretching, in agree- 
ment with experiments 0, 0, Il0| • Our explanation for 
Al contacts is in qualitative agreement with Ref. but 
differs from Ref. |l2( , where the displacement of the cen- 
tral atom from the z-axis is important. Equation J3J also 
predicts that G has an upper bound of Go and a lower 
bound of 0.75 G n, in reasonable agreement with the his- 
togram data [ill25|. and that T (sp contribution) in- 
creases during the stretching, while T\ and T_i decrease. 
Equation Q and its implications are the second main 
result of this paper. 

Pb contacts. - Next we apply this approach to Pb con- 
tacts, where 6p orbitals mediate the transport. In an iso- 
lated Pb atom, the spin-orbit coupling splits them into 
P1/2 an d P3/2 with the former lower in energy by about 
0.9 eV |2Q. In a Pb contact with the 27r/3 rotational 
symmetry, its channels then can be classified as (J, rrij) = 
(1/2, ±1/2), (3/2, ±3/2), and (3/2, ±1/2), where j is the 
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total angular momentum and rrij is its component along 
the z-axis. Since each pair of channels (j, irrij) are de- 
generate, the ± sign can be regarded as a fictitious spin 
variable for the channels. Then we have three transport 
channels, each with the "spin" degeneracy two, and ob- 



tain Ti 



2 J Ef dE \p( j>±rn .) if ,{E) - pQ i±m .^ (E)]. 

Since N^ s = 2 for a Pb atom (2 electrons in the 6s 
orbital are ignored), one obtains, instead of Eq. 
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sin 2 (7rA% | mj |)/2), where AN { 



U,\ m i\) 



{h\ m i\ 



\AN ( . 



= 2. 



(5) 



The role of the constraint (J5J becomes clear in the large 
d limit, where the orbital coupling to electrodes is weak 



and \AN, 



(3/2,3/2)1 
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tion (JSJ) then predicts T(3/2,3/2) — T( 



l A ^(i/2,i/2)l- Equa- 



(3/2,1/2) 



< T, 



(1/2, 



which agrees with a channel-resolved measurement 10] . 
Note that except for the difference in N' (or N^ s ), Pb 
and Al contacts are very similar. The difference in N' is 
however crucial; from Eq. Q, one obtains 



G_ 



r + cos ( ir AN, 



(3/2,3/2), 



(G) 



which differs from Eq. We next consider the behav- 
ior of AiY(3/ 2l 3/2) • Since the (3/2,3/2) channel is least 
favored by both orbital and spin-orbit couplings to elec- 
trodes, |A./V( 3 /2.3/2)| should be smaller than 2/3 [Eq. JSJl] 
and decay during the stretching. From Eq. JJjJl, G/Gq 
then decays monotonically during the stretching with the 
upper (lower) bound of 9/4 (0) , which is in reasonable 
agreement with experiments |10] . For realistic values of d, 
on the other hand, the orbital coupling will lift the degen- 
eracy \ANi 



(3/2,3/2) 



= \AN, 



(3/2,1/2)1 



However this effect 



turns out to be rather minor: A recent calculation [lfl ] 
suggests that already at d = 1.2do.pb, where do,pb is the 
bond length of bulk Pb, ^^(3/2,3/2)] and A7V (3/2 ,i/2) I 
(< 1) are significantly smaller than | AA^/2,1/2) (> !)• 
Then the predictions of Eq. JBJ should remain qualita- 
tively valid at least for d > 1.2e?o.pb- We also note that 
in the range of d, where \AN^/ 2 ,i/2) I > 1, n °t only G but 



also each individual Ti 



U,\ m j\) 



decreases during the stretch- 



ing due to Eq. JSJ. This part on Pb contacts constitutes 
the third main result of this paper. 

Resonance in Al contacts: Ab initio calculation. - 
While the resonance formation is already verified for 
Na [3E1E3 and Pb 10] contacts, it remains rather 
unclear for Al contacts |2j| ■ We thus perform first- 
principles calculation for an Al contact with the geometry 
in Fig. \l\(L = 3). Pseudopotentials and real-space 
multigrid method |29| within the local-density-functional 
approximation are used. Each electrode is modeled by a 
fee cluster of M (= 69) Al atoms. Small cluster calcula- 
tions are useful to test the resonance formation although 
the predicted resonance positions may be wrong due to 



2 o 

CO 

' 2 

1/3 

o 

9 1 

o 

2 



(&)d = d ,L = 3 



(b)d = 1.3 d , L = 3 



(c)d = d ,L = l 



(d)d = 1.3d , L = 1 
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FIG. 2: Projected DOSs (PDOS) onto the central atom of 
Al contacts with different d and L. Solid and dash-dotted 
(dotted and dashed) lines represent the PDOSs for the m — 
(rn — ±1) channel with even and odd parity, respectively. Ef 
is aligned at eV and the PDOS for m = with odd parity 
is magnified by factor 2 for clarity. 



the finite size effect (at best < 1 eV in our estimation) . A 
supercell geometry containing the whole system is used 
and electrodes in neighboring supercells are separated by 
more than 10.58 A, so that inter-supercell interactions 
are negligible. Using a grid spacing of 0.44 A, the total 
energy is converged to within 10~ 7 Ry. For the M = 69 
cluster, single particle level spacing Se is estimated at 
~ 0.07 eV near Ef- To extract properties of infinite 
metallic electrodes from the finite cluster calculation, the 
Fermi-Dirac level broadening Et should be larger than 
Se. Using Et = 0.4 eV (> 3Se), we obtain the projected 
DOSs (PDOS) of different m states onto a region around 
the central atom (a sphere with the radius of 0.56 do) 
and decompose them into even and odd parity compo- 
nents with respect to the inversion symmetry. Here, do 
(= 2.86 A) is the bond length of Al fee metals. 

The resonance assumption is tested for d = do and 
1.3 do (Fig-EJ- For d — 1.3 do, the stretching distance is 
2(tf — do) = 0-6 do, which is smaller than the experimental 
stretchingdistance of 0.7 do to 1.7 do before the contact 
breaking HUGH- Figure|2b) shows the PDOSs for d = 
1.3 do- In the m = channel, two well-defined resonances 
develop, one with even parity (almost fully filled) and 
the other with odd parity. In the m = ±1 channel, one 
odd parity resonance develops. Thus for d = 1.3 do, the 
resonance assumption holds for all channels. When d 
is reduced to d [Fig.|2[a)], channels respond differently 
due to the difference in the coupling strengths between 
the a and tt couplings. Whereas the resonance in the 
to = ±1 channel remains well-behaved, resonances in the 
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m = channel become broader and make nonnegligible 
overlap at Ep. Thus not all channels are in the clear 
resonant regime. Recalling that the resonance widths 
in Fig. [21 can be overestimated by Ex, we may safely 
conclude that Al contacts enter the clear resonant regime 
somewhere between do and 1.3 do- For do < d < 1.3do, 
the correction to G from the m = channel being not in 
the clear resonant regime is — sin(7r?7) sin[7r(ry + 2AiV±i)], 
where r\ = 2 — 2R m= o. e - Thus for r\ <C 1, the correction 
is not significant and the general trend of rising G during 
the stretching is not altered. 

The PDOSs are calculated for the L = 1 geometry 
model of contacts as well, where the two electrodes share 
a common apex atom. The PDOSs for d — do [Fig. GIc)] 
are essentially structureless, implying that the sharp tip 
geometry alone does not induce resonances at Al apex 
atoms, in contrast to the result for the Na (111) tip ge- 
ometry 0, ^| . But the resonance assumption holds for 
d = 1.3 do [Fig. Etd)] and a conclusion similar to that for 
the L = 3 geometry can be made. 

Discussion and conclusion.- We first note that al- 
though contacts with a symmetric shape are used for 
analysis, deviations from the symmetric geometry such 
as disorder in atomic positions do not affect the result 
qualitatively, as demonstrated explicitly in Refs. fTlt f]~3| . 
Secondly, our analysis can be used to study effects of L 
variation as well. For monovalent contacts, for example, 
the resonance analysis predicts different conductance be- 
havior for even and odd L due to the difference in N^ s . 
This prediction is consistent with the oscillation of G with 
L in Na contacts (first- prin ciples calculations (l4Tll5j |'). In 
Au, where L is known |3JJ| to vary duringthe stretching, 
a similar oscillation was indeed observed [a, 0| , although 
its amplitude is smaller than the predicted value for Na 
probably due to the stronger coupling in Au • Thirdly, 
our analysis can be also extended to contacts made of sd 
metal, which have five channels and thus could be more 
complex. It was recently reported |13 that a one-atom 
Pt contact has effectively three open channels (sd z 2, d yz , 
and d zx characters) at the Fermi level under zero bias. 
Here, one can find N' eS = 4 since 6 among 10 valence elec- 
trons of a Pt atom participate in the almost fully filled 
sd z 2, d xy , and d x i_ y i states Then, our resonance 

analysis leads to G/G ^ 9/4 - (1/2 + costtAA^/^ J 2 
and \AN sdz2 \+2\AN dyz/d J = 4. Because \AN dy ]/ d J is 
larger than 4/3 and increases (at best up to 2) during the 
stretching 33], G decreases during the stretching with 
the upper bound of ~ 9/4Go, in good agreement with 
Ref. [22|. Fourthly, we remark that our approach may 
be also useful for a wider class of systems, for example, 
symmetric molecules coupled to electrodes [34|. For this 
purpose, the charge neutrality condition (01 should be 
generalized to take into account possible electron trans- 
fer between molecules and electrodes arising from their 
electron affinity difference. 

To conclude, under the single assumption of resonant 



transport, it is demonstrated that qualitatively different 
conductance behavior in various stretched contacts can 
be explained in a unified way in terms of the number of 
electrons participating in partially filled resonance states, 
orbital shapes, and charge neutrality. 
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